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This is the fifth of a series of reports 
covering an investigation of the general 
instability problem by the California 
Institute of Technology. The first five 
reports of this series cover investiga- 
tions of the general instability problem 
under the loading conditions of pure bend- 
ing and were prepared under the sponsor- 
ship of the Civil Aeronautics Administra- 
tion. The succeeding reports of this se- 
ries cover the work done on other loading 
conditions under the sponsorship of the 
National Advisory Committee for Aeronautics. 



SUMMAEY 



This report summarizes the work that has been car- 
ried on in the experimental investigation of the problem 
of general instability of stiffened metal cylinders sub- 
jected to pure bending at the O.I+f, This part of the 
investigation included tests of 46 sheet-covered speci- 
mens. The most significant result was the determination 
of a new design parameter for the case ©f a stiffened 
metal cylinder subjected to pure bending. 

INTRODUCTION 

It is intended to give in this report a summary of 
the results of the experimental investigation of the 
general instability of stiffened metal cylinders sub- 
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jocted to pure bending. The earlier work at C.I.T. on the 
problem of general instability of stiffened metal cylin- 
ders has been reported in references 1, 2, 3, and 4. A 
total of 46 sheet- cover ed specimens has been tested. 

The investigation was planned so as to include most 
of the essential variables involved in the problem and 
was sufficient in scope to determine a suitable design 
parameter. Although this parameter was obtained for spec- 
imens subjected to pure bending, it is felt that it may 
be possible through suitable modifications to use the 
same parameter for combined loadings. The experimental 
data are also sufficient to allow a thorough check on 
theoretical methods for predicting general instability* 
A comparison of the predicted and experimental results is, 
therefore, included in this report. 

In addition to the above results, an account is also 
given of the experimental studies which have been made to 
give a better understanding of the mechanism involved in 
the buckling phenomena of stiffened cylinders. The body 
of this report, therefore, consists of four parts; namely, 

(a) Experimental investigations on the failure of 

metal cylinders 

(b) A comparison between theoretical predicted gen- 

eral instability stresses and the experi- 
mental results 

(c) A discussion of the mechanism involved in the 

buckling phenomena of stiffened cylinders 



(d) A new design parameter 



EXPERIMENTAL INVESTIGATION 



In starting this research program, it was felt that 
a complete study of ''the pure bending phenomena would form 
a desirable background for the more complicated problem 
of bending plus shear and for the more general combined 
loading conditions. With this idea in mind, the research 
to date has been confined to a systematic investigation 
of cylinders subjected to pure bending. 



To obtain a better understanding of the amount of 
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experimental research, work which can new be considered as 
Completed, it might he well to hriefly review the varia- 
bles involved in the prohlem of general instability. The 
variables involved may be divided into two classes - those 
dealing with the geometry of the structure and those in- 
volving the sectional properties of the stiffening ele- 
ments as well as the sheet covering. 

The geometrical variables involve: 

(a) a variation of longitudinal spacing 

(b) a variation of frame spacing 

(c) .a variation in diameter 

(d) a variation in length 

• The second group of variables involves: 

(a) a change in the section properties of the lon- 

gitudinals 

(b) a change in the section properties of the frames 

(c) a change in the thickness of the sheet covering 

^The first group of variables has been thoroughly in- 
vestigated inasmuch as a large number of specimens having 
2.o3-, b. 06-, and 10.12-inch longitudinal spacings, and 
1 ~' 2 ~ • 4 7> 8 ~i a nd 15-inch frame spacings have been 
tested. Although. the second group of variables is not yet 
completed, a sufficient number of the variables have been 
investigated to form a basis for the analysis of general, 
instability. Three different frame sizes have been in- 
vestigated; that is, frames 5\, F 5 . and f e (fig. i), and 
also two different sheet thicknesses, 0.010- and 0.016- 
inch have been used. Six specimens, having a 20-inch diam- 
eter have been tested; all other specimens are 32 inches 
m diameter. Of the former', two had an L/D ratio of 2.6, 
while in all other cases the L/D ratio was 2.0. The lon- 
gitudinals have not been varied. However, by varying the 
sheet thickness, a certain variation in the section pron- 
e ^ le f. 0f tbe lon Bi'tudinals is obtained, since a certain 
effective width of sheet which varies with the sheet thick- 
ness is assumed to act with each longitudinal. 

These variables which have been investigated have been 
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found sufficient to establish a suitable design parameter 
for structures subjected to pure bending. Hence, if we 
confine ourselves to the pure-bending phase of the inves- 
tigation, it is felt that future work will only involve a 
number of check specimens in which the longitudinal sec- 
tion properties are varied. Also, an additional number 
of specimens should be tested to further investigate the 
influence of local instability of the frames. A few ad- 
ditional specimens will be tested to further check the 
length effect . 

A complete list containing all the necessary data of 
the sheet-covered specimens which have been tested to date 
is given in table I. It should be noted that the maximum 
unit compressive deformation Cj/B 1 at failure is tabulat- 
ed rather than the compressive stress, where G is the 
compressive stress and B 1 is the effective modulus - the 
reason being that in many cases the compressive stress at 
buckling exceeds the yield point of the material, and the 
difficulty involved in determining the basic stress in 
such cases makes it desirable to use the unit strain 
rather than stress. The nondimen si onal quantity cr/E 1 is 
also more desirable from the standpoint of determining a 
nondimensi onal design parameter. 

The method of measuring the unit strain at failure 
was discussed in reference 4 and needs no further dis- 
cussion. It might be of interest, however, to point out 
that the accuracy of the method is somewhat dependent on 
the type of failure. Until local bending occurs the 
method is auite accurate; but if bowing between frames 
takes place, then the test data must be corrected for bend- 
ing and, consequently, are subject to the inaccuracies of 
such calculations. General instability failure usually 
occurs by a sudden collapse of the structure, whereas in 
panel instability the failure in most cases is gradual; 
that is, the applied bending moment will increase after 
buckling of the longitudinals between frames occurs. Also, 
as discussed in reference 4, the sheet will reinforce the 
longitudinals after the amplitude of the longitudinals 
becomes sufficiently large - thereby causing a higher and 
less clearly defined ultimate load. It is therefore much 
more difficult to determine the ultimate load as well as 
the longitudinal deformation of specimens which fail by 
panel instability. 

Figures 2 to 26 indicate the unit longitudinal defor- 
mation as a function of the applied bending moment for 
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the sheet-covered specimens. All specimens failed either 
by panel or general instability with the exception of 
specimen 62, which failed in tension. Failure of this 
specimen was precipitated by a bearing failure of the sheet 
at the joint, thus causing the entire tension load to be 
carried by the longitudinals. 

It was pointed out in reference 4 that no correlation 
could he obtained between the sheet-covered specimens and 
the wire-braced specimens previously tested. Hence, data 
pertaining to the wire-braced specimens have not been in- 
cluded in this report. 

A number of photographs (figs. 43 to 48) showing the. 
various types of failure are included. 



TKSOSST I CALL Y PREDICTED GE1TSRAL INSTABILITY STRESS 



A detailed discussion of the analytical methods which 
are available for calculating the general instability 
stress of stiffened cylinders is given in reference 1. As 
pointed out, the available experimental data were too mea- 
ger to indicate whether any of the methods could be satis- 
factorily used for the prediction of general instability. 
Sufficient experimental data are now available to deter- 
mine the merits of each of the methods by comparing the 
predicted general instability stress with the experimental 
values. Each method will be discussed separately. 



m 

t i on 



Hof f 1 s Method 

The general instability stress is given by the equa- 



A x (Ad) 
where 



TT 2 S v I 

(1) 



E x Young's modulus of elasticity, pounds per square inch 

I x moment of inertia of a longitudinal plus an effective 
width of sheet (See section on a Hfew Design Param- 
eter, for effective width calculations.) 
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A x cross-sectional area of. a longitudinal 
d frame spacing 

The nondimensi cnal quantity X is a function of the 
structural coefficient A and the number of frames m 
involved in the failure. The relation between A and A 
for constant values of ra is shown in reference 2. It is 
seen that, for a constant value of A, the value of A 
increases as the number of frames involved per wave length 
increases. This increase in A produces, according to 
equation (l), a decrease in the general instability stress. 
In applying Hoff 1 s method to a practical design problem, 
it would be necessary to consider the largest number of 
frames which could possibly fail, hence give the lowest 
possible stress. However, in making the comparison be- 
tween the unit deformation a/E 1 , as calculated by equa- 
tion (l), and the test value, the actual number of frames 
which failed was used. A comparison of the experimental 
and calculated values for a large number of specimens is 
shown in table II. The calculated values range from 11 to 
83 percent of the experimental values. If the total num- 
ber of frames in the test specimen were to be used in the 
calculations, the calculated values would be lower. It may 
therefore be concluded that, although in all cases Hoff s 
method does predict general instability, the method is too 
conservative for design purposes. 



Dschou 1 s Method 

Dschou obtained for the minimum general instability 
stress an equation of the form 




which can be written as 
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wh ere 

p x radius of gyration of a longitudinal together with 
the corresponding portion of the sheet covering 

py same quantity for the circumferential frames 

t x = t + ~ 2l s equivalent shell thickness in the axial 

direction 

A y 

ty « t + — , equivalent shell thickness in the circum- 
ferential direction 

t sheet thickness 

A x area of a longitudinal 

ky frame area 

The same equatien is reproduced in T imc shonko 1 s theory of 
elastic stability (reference 6). 

It is not definitely stated by either Dschou or Timo- 
shenko whether, in calculating I x , I y or, similarly, 

pt pi the total width of sheet between longitudinals 
~ - $ 

and frames should "be used or an effective width of sheet. 
Calculations have been made, using first the total width 
of sheet "between longitudinals nnd frames and then using 
an equal effective width of sheet with the longitudinals 
and frames. The effective width calculations are dis- 
cussed Later in this report. The results of these calcu- 
lations are shown in figures 27 ^.nd 28, The calculated 
values in all cases are considerably higher than the ex- 
perimental values. It is also evident that the expression 



2 
R 



is net suitable as a parameter, since the scatter of the 
experimental points is much larger than can be contrib- 
uted- to experimental inaccuracies. 



/ p * P * !* 



X 



y 



J 



P x -+P. 



2 *X 



(4.) 
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Kyder 1 s Method 

The application of this method requires the calcula- 
tion of three parameters, K 2 , K 3 , and K 4 , depending on 
the geometry of the structure and on the section proper- 
ties.' Curves of X 2 and K 3 are not available for val- 
ues of K 4 < 0.20; hence, in all cases where K 4 was 
less than 0.19, it was felt that the curves in which 
K 4 ss 0.20 would no longer he valid, and these specimens 

were therefore not considered. The results of the calcu- 
lations are shown in table III, where Pjj^OB are tiie 
calculated general instability stresses. The ratios of 
the calculated to the experimental values are also shown. 
The experimental stress in this case is not a true stress 
hut an equivalent stress which is obtained from the meas- 
ured unit deformation at failure under the assumption 
that the modulus of elasticity remains constant beyond 
the proportional limit of the material. The results indi- 
cate, with the exception of three specimens, that the cal- 
culated values are all considerably higher than the ex- 
perimental values. It should bo pointed out, however, 
that in practically all cases the value of K s is less 
than 0.2, which makes it extremely difficult to obtain 
the value of t e /p s from the curves since in this range 
the value of t e /p s is very sensitive to a small change 
in K a . 

Taylor 1 s Method 

The general instability stress was calculated f6r a 
number of specimens and is shown in table IV. The calcu- 
lated stresses are of the order of five to twelve times 
higher than the experimental values. The given calculated 
values are from Taylor l s approximate equation. It was 
found that the exact equation gave values only slightly 
different (ar>prox. 2 percent) from that of the approxi- 
mate equation, and the simplicity of the latter makes it 
preferable for calculations. 

A STUDY OF BUCKLIUG PHENOMENA OF STIFFENED CYLINDERS 

In the course of the experimental investigation, it 
was observed that in all cases of general instability the 
last specimens showed a definite preference for failing 
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radially inward. Also in measuring the "normal restraint 
coefficient," it was observed that the radial deflection 
was a nonlinear function of the applied load. (See fig. 
32.) Furthermore, if the load is applied radially out- 
ward, the slope of the 1 oad- de f 1 e ct i o n curve increases, 
while for a load applied radially inward the slope de- 
creases. This simply means that the structure has greater 
rigidity in the outward direction than in the inward di- 
rection. 

Considering the stiffened cylinder as a whole, it is 
immediately evident that it is anisotropic and the influ- 
ence of one member upon the other is extremely difficult 
to determine. The most elementary concept of the problem 
would he that of a column supported "by continuous and con- 
centrated elastic supports - a longitudinal stiffener be- 
ing the column, the sheet covering providing the continu- 
ous elastic support, and the frames the concentrated elas- 
tic supp or t . 

The problem of a column elastically supported has bee 
discussed by H. Zimmerman for the case of concentrated sup 
ports, and by F. Sngesser and others for the case of a dis 
tributed support. In all cases the investigations have 
been confined to elastic supports exhibiting a linear 
force-deflection relation. 

Since the frames and the sheet have the characteris- 
tics of a nonlinear elastic support, it was felt that it 
would be of general interest to investigate the effect on 
the load-deflection relation of a column supported by a 
nonlinear elastic element. Inasmuch as rings or frames 
are known to have the desired nonlinear properties, a thin 
semicircular steel ring, as shown in figure 43, was used 
as the elastic support. For simplicity, columns having a 
simple, rather than a distributed, support were tested. 
It may be of some interest to consider first the elastic 
behavior of a semicircular ring. Designating the radial 
load by P, and the corresponding radial deflection by 6, 
the curves of figure 29 indicate that if the load is ap- 
plied radially inward, the value of P/S, the elastic 
constant, decreases with increasing deflection. When the 
load is directed radially outward, the value of P/S 'in- 
creases with increasing deflection. Obviously, then, if 
an initially straight column, is supported by such an ele- 
ment, or elements, it may be expected that it would have 
a tendency to buckle in the direction of "decreasing P/8 
or, if due to an initial deformation in the direction of 
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increasing P/fi the buckling starts in that direction, 
then at some deflection a sudden ".jump" may occur in the 
direction of decreasing P/S. . 

Tests #ere conducted on columns 0..090 inch thick by 
0.375 inch vide and 19 inches long. These columns were 
cut from 24S-RT alclad sheet stock. The steel ring was in 
all cases 8 inches in diameter and 0.008 inch thick. The 
test apparatus and method of testing are illustrated in 
figure 43. 

.The. results of those tests are indicated in figure 
30, where the ratio of the column load to the Suler load 
is plotted as a function of the ratio of the normal de- 
flection 8 at the center of the column to the column 
length I . 

Considering the results of figure 30, it is seen that 
the elastic support increases the buckling load of the 
u straight" column (upper curve) to nearly 3.5 times the 
Suler load. This load is reached at a relatively small de- 
flection. As the deflection 6 increases, the decrease 
in load is at first quite rapid, then more gradual as the 
deflection "becomes larger, and may approach a minimum at 
large deflections. It was not possible to reach very 
large deflections "because of plastic failure of both the 
rings and the columns. The lower curves in the same fig- 
ure indicate the effe.ct of initial deflections in which 
the column was rolled to approximately the form of a. 
half- sine-wave t the maximum initial deflection being des- 
ignated .by 8 0 . These curves show that, .with increasing 
initial deflection, the maximum load decreases and occurs 
at increasingly larger deflections. In all cases the load 
sustained by the column tends to approach at. large deflec- 
tions the "minimum load" of the straight column. Thus, 
in case of a nonlinear support, any initial imperfections 
of the specimen will appreciably lower it.s buckling load. 
To illustrate the contrast between the column with a non- 
linear elastic support and one with a linear elastic sup- 
port (a coiled spring), a number of tests were conducted 
on the same type of column. As may be seen from the 
curves of figure 31, the columns with an initial deflec- 
tion in all cases approach the maximum load of the straight 
column. 

The most striking features of the column supported by 
a nonlinear elastic element are: first, as the deflection 
increases the load decreases; secondly, there are two or 
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more possible configurations of equilibrium for the same 
load; one corresponds to 6 = 0 t . the others to 8 > 0. 

Since the longitudinal members in a stiffened cylin- 
der will "behave in a manner similar to the column with a 
nonlinear elastic support (fig. 32), the buckling charac- 
teristics of a stiffened cylinder must be similar to those 
of the nonlinearly supported column. Therefore, any the- 
ory which is based on the assumption of small deflections 
will probably give a much higher load than that actually 
observed, since In a test specimen initial imperfections 
will prevent the specimen from reaching the maximum load, 
and will cause it to fail at some lower load depending on 
the degree of imperfection present in the specimen. This 
expectation is verified by the comparison between the ex- 
perimentally observed stresses and the predicted stresses 
as obtained from the theories developed by Taylor, Dschou, 
and Ryder, which are essentially based on small deflec- 
tions and linear restraints. 

A theory based on large deflections would be extreme- 
ly difficult because of the nonlinear load-deflection re- 
lationship discussed in this section. This difficulty 
makes it desirable to attack the problem by experimental 
methods which would lead to the development of a suitable 
design parameter. 

A JJ2W DESIG-1T PARAMETER 



The buckling phenomena of stiffened cylinders and the 
attending difficulties of a theoretical solution have been 
discussed in the previous section. It was hoped, however, 
thcat it would be possible to obtain a design parameter 
which would determine the critical value of cr/3 1 with a 
satisfactory accuracy for design purposes. 

With this thought in mind, a number of specimens in 
which only the geometry of the structure was varied were 
first tested. The results of these tests are shown in 
figure 33, where the unit strain cr/3 1 is plotted as a 
function of the longitudinal spacing b for constant 
values of d/b - referred to as the reinforcement aspect 
ratio,, The curve corresponding to an aspect ratio of 
1.58 was considered as the master curve and the remain- 
ing curves nre an expansion of" this curve; hence they fit 
best by changing the horizontal scale. 
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A log-log plot, of the expansion factor as a function 
of d/b indicated that this factor varied approximately 

as /d/b. In-fact, plotting the unit strain a/E r as a 
v ' 

function of ,^/bd , a single curve, as indicated in fig- s 
ure 34, is obtained. A log-log plot of the unit strain 

cr/E 1 as a function of J. bd showed that E 1 /(J varied 
approximately as -/bd. This, then, gives the influence 
of the geometrical variables on the general instability 
f ailur e . 

The next question is: In what manner does the radius 
B and the section parameters p and p enter into the 

V 

design parameter sought? By analogy - with the buckling 
of unstiffened cylinders - it can he expected that, for 
identical values of b, d, p x , and p the reciprocal of 

the critical value of the unit strain varies linearly with 
E. This surmise was checked by testing specimens with two 
different radii (H = 1? in. and 10 in.) and was found to 
be correct. Hence, it vas concluded from these results 
that the design parameter has the form 

4 

V bd R 



f (p , P ) 
x r y 



on 
a 



From dimensional reasoning, it follows that the functi 
f (p x ,p v ) must have the dimension of the 3/2 power of 

length. The simplest assumption for the function which - 
determines the influence of the section parameters, p x 

and p , is that it depends on the geometrical mean value 

y 

J p p only. It was found that, in the range covered "by 
x y 

the experiments and by the use of reasonable methods for 
the calculation of the effective width of the sheet, .this 
assumption appears justified. Thus the design parameter 
appears in the form 



v 



hi 7 
or 



p ) J P x K J p x p v 

x y v x y x y 
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Concerning the width of sheet to "be used with the Ion 
gitudinalfl and frames, Hof f (reference 7) and Cox (refer- 
ence 5) have pointed out that the effective width for sta- 
bility is not the same as that based on the load-carrying 
ability of the sheet. In fact, Cox states that the effec- 
tive width associated with buckling phenomena is propor- 
tional to the rate of increase of the apparent stress to 
the rate of increase of the actual stress, rather than to 
the chord value as illustrated in figure 3 (reference 5). 
Assuming that Harquerre 1 s effective width equation 

b 

is sufficiently accurate for our purpose, then the appar- 
ent stress G" is given "by the equation : ., 

a „ st st J a st 

wher e 

longitudinal stress, pounds per square inch 
(~ edge stress) 

0"- ."buckling stress of 'the sheet, pounds per square inch 
c 

The effective width for stability is then 



3v e _ dg a _ 2 / u c 

* " 3 y^t • 

In all calculations involving the effective width of sheet 
acting with the longitudinals, this value has been used 
with the exception of the calculations involving Ryder's 
work, since in his work the effective width is definitely 
speci f i ed . 

The amount of sheet acting with the frames is diffi- 
cult to evaluate by analytical means; trial calculations 
indicated that the "best results were obtained when the to- 
tal '< r idth of sheet between frames was used. For this rea- 
son, py ■ was calculated assuming the entire width of 
sheet to be effective. The variation of X x , I y , p x , 

and p a's a function of the effective width is shown in 
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figures 38 to 41, where I, 



and 



Po 



is the moment of in- 



ertia and the radius of gyration of the member alone. 

In some cases a further correction, which is discussed 
in the appendix, is to he applied to the radius of gyration. 



The unit strain a/E ' 



plotted as a function of 



at failure (general instability), 

A i 

S 



/- 

J P 



hd 



i s shown in 



figure 35; this curve 
purposes as the slope 
high values of or/E 1 . 



VP p 

y * 

is not very convenient for design 
of the curve "becomes too steep for 
A more desirable presentation is 



to plot cr/E 1 as a function of 



y 



/ 



Pv 



R 



J 



bd 



The 



results are indicated in figure 36. As may be seen, the 
experimental values lie close to a straight line with the 
exception of three points. Two of these points corre- 
spond to tests on the channel section frame. As these 
frames failed by local instability of the outstanding legs 
(see fig. 44), it may be expected that the results would 
be low for the channel does not develop the strength corre- 
sponding to the calculated value of p . This' can be il- 
lustrated by considering the behavior of an open-section 
column subjected to an axial thrust* The critical column 
stress is given by the equation 



TT 



If the critical column stress is higher than the local in- 
stability stress, the column will fail at the lower stress 
This essentially means that an effective value of p can 
be obtained corresponding to this lower stress, that is, 



Peff _ 



(7 



J 



f ai lur e 



It seems reasonable to assume that the frame behaves 
in a similar manner and that the effective radius of gy- 
ration is lower, due to the local instability failure. 



This new design parameter, 



I 



j~? p y p p 

v x y x y 



bd 



i is the 
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ratio between the radius and a quantity having the dimen- 



sions of a length defined by J pgp^ J • ; y . Hence, it 

appears that this length is proportional, to some extent, 
to an equivalent thickness of the reinforced cylindrical 
structure . 

The failing stress of an unstiffened cylindrical shell 
of radius R is given "by 

k i ( 5 ) 

E 1 B 

where t is the thickness of the shell and k is an em- 
pirically determined constant. Replacing the thickness 
t by the radius of gyration p of a strip of the shell 
of unit width, equation (5) can he written in the form 



cr 

TT; « 



= k (e) 



For purposes of comparison, the failing stress of the 
reinforced cylinders investigated can "be expressed as a 
function of the new parameter in the form 

* . k , ^ /EE (7) 

E 1 R 

where k 1 is a numerical constant. Figure 36 shows that 
equation (7) covers, with a fair approximation, the cases 
in which general instability occurred. 



By introducing tho geometrical mean value J p x p v , 
it can "be assumed that the influence of the anisotropy of 
the structure is approximately taken into account. Then, 
"by comparison of equations (6) and (7), an effective radi- 
us of gyration can he defined as 



si y^fx iW = „ /7"p~ (s) 

k J bd r x y ^ * y 



where cp = — / £ can "be considered as a correction 

k V "bd 
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V>x Py a 

factor to the ratio 5 • Now — is the slenderness 

R Px 

ratio of the longitudinal, taken as a column between two 
frames and, similarly, — is the slenderness ratio of the 

Py 

frame taken as a column between two longitudinals. If the 

geometrical mean value of these two slenderness ratios is 
defined by A, that is, 



A = /-i4_ (9) 



then the correction factor can "be expressed as 

do) 



kj 1 



It is evident that if the stiffened shell were to "be con- 
sidered as an equivalent shell, in which all the materi- 
als are uniformly distributed, then the only appropriate 



/ P P 

parameter to he used is . The appropriate param- 

B 

eters which enter into the problem of 'buckling of a truss 
are the slenderness ratios which appear in the quantity 



X zz / — — . The fact that the empirically derived rela- 

/p p~~ 

tion given here involves "both parameters, " and A, 

indicates that a stiffened cylinder cannot he treated ei- 
ther as an equivalent cylinder of uniform thickness or as 
a cylindrical truss. . 'Therefore, the work of Dschou, 
Taylor, and Timosh-enko, which is based upon the concept of 
an equivalent uniform shell, cannot be expected to give 
correctly the general instability stress of stiffened cyl- 
inders* The error is particularly large if the flexural 
rigidity is concentrated in frames of large cross-sectional 
ar ea . 

It is of interest to calculate the numerical value of 
the factor rp , using for k 1 the values obtained from 
the experiments of this report, and for k the value 
k ss 0.3, which is a reasonable average value for unstiff- 
ened shells in the range of (j/B involved. The values of 
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4 



J 



1 Fx p y 

bd " 



vary from 0.0598 to 0.219, as shown in table V; 



the average slenderncss ratio A varies from 280 to 21* 
The values of cp corresponding to these limiting cases 
are cp = 0.2154 and op = 0.932, The lower limit is ob- 
tained for a structure (specimen 34) with a 10. 12-inch 
longitudinal, a 16-inch frame spacing, and an average 
slenderness ratio of 280. The upper limit refers to spec 
imen 55 xvith a 2.53-inch longitudinal, a 2-inch frame 
spacing, and an average slenderness ratio of 21. 

It is seen that for structures in which the stiffen- 
ing elements are widely spaced, it is necessary that a 
multiplying factor much smaller than unity he applied to 

the quantity J p v P„ • In the case of close spacings the 
^ x y 

factor is of the order of 1. 

It is realized that this discussion does not repre- 
sent a complete theory of the "buckling of reinforced cyl- 
indrical shells subjected to pure "bending. In particu- 
lar, it cannot "be expected that the use of geometrical 
mean values will take care of the anisotropy of the struc 
ture with sufficient accuracy in all extreme cases; for 
example, if the cylinder is stiffened in only one direc- 
tion. It is hoped, however, that the relation developed 
in this report will cover most of the structures occur- 
ring in practical design. 

•The practical application of the design parameter is 
relatively simple . Since the designer will have .the nec- 
essary data to compute the numerical value of the parame- 
ter, the general instability stress can he immediately oh 
tained from the curve of figure 37. In view of the fact 
that general instability causes a complete collapse of 
the structure, it is recommended that allowance he made 
for ample margins, and in no case should the allowable 
stress exceed the boundary line indicated in figure 37. 
As previously stated, when the frames fail by local insta- 
bility the general instability stress is considerably 
lover than that defined by the curve of figure 37. Hence 
until further data are obtained on this type of failure, 
the curve of figure 37 should be used with caution when 
local instability of the frames i s likely to occur. 

If a reinforced cylinder fails by panel instability, 
the buckling stress should in all cases be lower than the 
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stress necessary to cause a general instability failure. 
This is obvious since, in order that panel instability may 
occur, it is necessary that the frame he sufficiently rig- 
id to maintain closely the shape of the structure at the 
frame. If the frame is not sufficiently rigid it will 
fail before the panel instability stress is reached Rnd 
result in a general instability type of failure. A number 
of the test specimens failed by panel instability. The 
unit strain at buckling is shown in figure 37 and, as 
pointed out, these values should and do lie below the 
curve defining general instability. 



Guggenheim Aeronautical Laboratory, 

California Institute of Technology, 
Pasadena, Calif., July 1940. 



COERSCTIOH FOR T KS EFFECT OF THE TORSIONAL RIGIDITY 
OF FRAMES 021 TEE STRENGTH OF CYLINDERS 



Owing to the particular method of attaching the stiff- 
oners to the frames, it is felt that the torsional rigid- 
ity of the frames acts as an additional restraint in the 
bending of the stiffeners. This effect will be apprecia- 
ble if a frame such as F 6 , which has a large torsional 
rigidity, is used. It is therefore necessary to correct 
for this effective increase in rigidity. A physical rea- 
soning immediately leads to the following expression for 
the effective moment of inertia of the stiffeners: 



APPEUDIX 



eff 



4- k C 



(1) 



wher e 




4 

effective moment of inertia of the stiffener, in. 



eff 



moment of inertia of the stiffener. "with effective 
width of sheet, in. 4 



C 



torsional rigidity of the frame divided by the 
shear modulus of the material, in. 
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k a nondimensional factor which gives the effective- 

ness of the fram.es in supporting the stiffeners 
in "bending 

The following analysis is an attempt to estimate 
this nondimensional factor k. Consider a wave of the 
Duckled shell, as shown in figure a. 



Jl 



- n 



n 



r 



T 

d 



I 

H 



i 



\ 



jLL 



MA 



I 



i ! ' ll \ 



-H 



Fran as 



\ Stiffeners 



?igure a. 



At the boundary of the wave the frames are not twisted 
"by^the "bending of the stiffeners. At the joint with the 
stiffener the frames are rotated an angular amount eaual 
to the slope of the deflection curve of the stiffener. 
If w is the deflection of the middle stiffener, then 
this angle is ~E at the f j oint v , ith thig stiffener. As 

a rough estimate, let the rate *£ rotation of the frame be 
~/l where I i s the average buckled wave length of the 
frame s • 
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The torsional rigidity of the frame is CG where G 

is the shear modulus, or 0 = , ^ r • *n order to 

2(1 + M») 

simplify the calculation, this amount of torsional rigid- 
ity is assumed to be distributed within the frame spacing 
d. Thus the distributed torsional rigidity is ^ ( x + ^ ) "a ' 

This distributed torsional rigidity multiplied by the rate 
of rotation of the frames gives the distributed moment 
acting on the stiffener due to the torsional rigidity of 
the frames. The total moment at a certain section due to 
this effect will be the integral of the distributed mo- 
ment from one end of the stiffener to the section con- 
cerned; or 

x 

P S C iy 1 . Y 

^ B I — — CLi. 

'torsion J g(l +-M») d dx I 



o 

E C w 



2(1 + |x) d I 

Thus the equilibrium equation for the stiffener is 



(2) 



SL dx^ 



C 1 



2(1 + MO d I 



r ± + P 



W =0 (3) 



Comparing this equation with the equation for Buler column 
with effective moment of inertia I Q ^ 



we have 



'ef f 

SI dfx + P W = 0 (4) 

So eff dx* 



1 0 L « El 

T 2 



EI st „ EI st 2(1 + ^ X st 

C I I 

where I is the wave length of the buckle; or 
rr 2 TT S I st eff ^ 1 CI 



L 2 I st 2(1 + |0 I 8t dl 



that is, 
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st 



eff _ l 



st 



2(1 + m,)tt S dl I 



st 



therefore 



st 



eff 



= I 



st 



LStt (l + \t) dl- 



Comparing this equation with equation (l), the value of 
the interaction factor k can he retained as 



k = 



2n s (l + hO H 



(6) 



using the values of I and L from experimental 
data, it is possihle to estimate the value of this factor 
which is found to oe of the order of 1. For the majority 
of the frames, the value of 
with the exception cf frame 



'■ s t is much larger than C 
Hence, the correction was 



only applied to the test specimens in which frame 



used. 



wa s 
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TABLE I 

Pure Bending Tests of Longitudinal - Frame Combinations 

All Longitudinals Si 



Test 


Frame 
No. 


Sheet 
Thick, 


Long. 
Spacing 
b 


Frame 
Spacing 
d 


Radius 


Unit 
Strain 




/O 
fx 






kr 


Type of Failure 


4./ ' 








• 010 


2. 53 


8 


15.92 


.00170 


2.120 


.1168 


.02335 


.01145 


2822 


General Instability 


Zo 






2.53 


4. 




. 00208 


1.783 


.1165 


.02677 


.01320 


2150 


n n 


27 






2. 53 


2 




.00288 


1.500 


.1165 


,02836 


.01380 


1730 


• 

« it 


28 






2.53 


16 




.00130 


2.520 


.1169 


.01875 


.01010 


3970 


Started by Panel Instability 
TinAl Failure General Instability 


Of* 






5.06 


16 




.00092 


2.998 


.1160 


.01875 


.01023 


4660 


Panel Instability 


30 






5.08 


8 




.00140 


2.520 


.1165 


.02335 


.01190 


3370 


General Instability 


OA 






5.06 


4 




.00190 


2.120 


.1169 


.02677 


.01325 


2550 


n « 








5.06 


2 




,00256 


1.783 


.1167 


.02836 


.01380 


2060 


n u 


34 






10.12 


16 




.00088 


3.560 


.1110 


.01875 


.00975 


5820 


Started by Panel Instability- 
Final Failure General Instability 


35 






10.12 


8 




.00120 


2.998 


.1120 


.02335 


.01160 


4120 


General Instability 


36 






10.12 


4 




.00164 


2.520 


.1122 


.02677 


. 01280 


3135 


n n 


37 






10.12 


2 




.00200 


2.120 


.1117 


.02836 


.01330 


2540 


n n 


38 






5.06 


1 




.00283 


1.500. 


.1164 


c 02801 


.01360 


1755 


ti ii 


39^ 




.010 


2.53 


1 




.00326 


1,261 


.1170 


.02801 


.01360 


1478 





TABLE I (COKTD»J. 



CO 



Test 


Frame 
No. 


Sheet 
Thick* 


Long* 
Spacing 
b 


Frame 

S racing 
4 


Radius 


Max. 
Unit 
Strain 






<°y 

7 






Type of Failure 






vW 


40 




.015 


2. 53 


16 


15.92 


, 00104 


2U520 


.1175 


.01723 


.00955 


4210 


Panel Instability 


41 






2.53 


8 




.00182 


2.120 


.1180 


. 02176 


.01140 


2960 


General Instability 


42 






2.53 


4 




.00215 


1.783 


.1184 


.02608 


,01310 


2168 


t? ti 


43 






2.53 


2 




.00280 


1. 500 


.1184 


.02879 


. 01410 


1693 




44 






5.06 


16 




.00066 


2. 998 


.1140 


.01723 


. 0093 5 


5110 


Panel Instability 


45 






5.06 


8 




. 00146 


2.520 


.1150 


. 02176 


.01120 


3580 


General Instability 


46 






5.06 


4 




.00167 


2.120 


.1151 


. 02608 


. 0128 5 


2630 


n it 


47 






5.06 


2 




. 00249 


1.733 


.1151 


.02879 


. 01380 


2060 


n n 


48 






10.12 


16 * 




.00048 


3.560 


.104 


.01723 


.00870 


6520 


Jrmei lnsxaDiiiuy 


49 






10*12 


8 




.00128 


2.998 


.1049 


.02176 


.01050 


4540 


General Instability 


50 






10.12 


4 




.00150 


2.520 


.1060 


.02608 


.01195 


3360 


w n 


51 


F 5 


• 015 


10.12 


2 


15. 92 


.00250 


2.120 


.1061 


.02879 


.01300 


2600 




52 




.010 


2.53 


16 


16.00 


.00180 


2.520 


.1168 


0.1288 


.0430 


938 


Panel Instability 


53 






2.53 


8 




.00216 


2.120 


.1165 


0.0811 


.0303 


1120 


n it 


54 






2.53 


4 




.003732 


1.733 


.1168 


0.0982 


c 03 500 


812 


General Instability 



TABLE I (CONTD. ) 



Test 
* 


Frame 
No, 


Sheet 
Thick, 


Long. 

Spacing 
h 


Frame 
Spacing 
d 


Radius 


Max. 

Unit 
Strain 


r ° 




9 






TVr>e of Fkilm»« 




mr h 


55 


Pi 


.010 


2.55 


2 


16 


.004775 


1.500 


.1168 


.11000 


.0378 


635 


General Instability 


56 






5.06 


16 




.00090 


2. 998 


.11.50 


.06330 


.0250 


1920 


Panel Instability 


57 








8 




.00228 


2.520 


.1150 


.08110 


.0300 


1345 


Ranel Instability 


58 








4 




.004350 


2.120 


.1160 


.09820 


.0350 


965 


General Instability 


59 


h 






2 




,006154 


1.783 


.1165 


.11000 


,0378 


742 




60 








8 




.00243 


2.520 


.1270 


.06340 


.0251 


1715 


ffenel Instability 


61 








4 




.004927 


2.120 


.1270 


.06400 


.0254 


1245 


General Instability 


62 






5.06 


2 






1.783 




> 06180 






Tension Failure 


63 






2.55 


4 




•004045 


1,783 


• 1278 


.06400 


.0255 


1045 


General Instability 


64 


P 6 




2.53 


2 


16 


.005100 


1^500 


,1278 


.06180 


.0247 


905 


N u 


65 


P 5 




2.61 


4 


10 


.003150 


1*798 


.1163 


, 02677 


.0132 


1362 


n ti 


DO 






£.61 


2 




. 004200 


1. 510 


.1168 


.02036 


.0138 


1095 


n it 


67 






5.22 


4 




.002950 


2.175 


.116 5 


.02677 


.0132 


16 50 


if ff 


68 






5.22 


2 




,003750 


1.798 


.1166 


,02836 


.0138 


1300 


t? n 


73 






2.61 


4 




.003200 


1.798 


.1162 


c 02677 


.0132 


1365 


ti n 


74 


*5 


.010 


2.61 


2 


10 


. 004400 


1.510 


.1164 


fr 02836 


.0138 


1095 


ff it 
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TABLE IT 

Calculation of the General Instability Stress 
by Hoff's Method 



Test N-o. 




• r.. 




A 


m 


LogA 


X 


% 




30 


8.46 


.0006676 


.000037 


447 


4 


2.65 


2.59 


.000561 


0.406 


31 


9*75 


.000668 


.00003624 


5,610 


8 


3. 558 


4.38 


.000666 


0. 3456 


32 


9.16 


.000664 


.00003655 


28 , 820 


15 


4.46 


7.85 


.000821 


0.3206 


35 


9.2 


.000815 


.0000464 


218 


5 


2.339 


2.8 


.000496 


0.415 


56 


9.65 


.000809 


.00004625 


1,740 * 


6 


3.241 


3.36 


.001365 


0.832 


37 


8.04 


.000823 


.00004632 


14 , 120 


16 


4.15 


8.5 


.000912 


0.454 


41 


4.595 


.000662 


.00003978 


826 


7 


2.917 


3.8 


.000219 


0.175 


42 


6.155 


.000643 


.0000370 


6,910 


11 


3. 84 


5.8 


.000365 


0.2116 


43 


6.5 


.0006395 
i — 


.0000387 


55„460 


16 


4.744 


8.4 


.000692 


0.28 


46 


5.84 


.000795 


.0000481 


410 


7 


2.613 


3. 72 


.0002742 


0.188 


46 


6.30 


.00079 


.0000478 


3,299 


11 


3.518 


5.74 


.0004575 


0.274 


47 


6.23 


.00079 


.0000478 


26,210 


16 


4.419 


8. 35 


.00087 


0. 358 


49 


6.775 


.000938 


.0000586 


198 


7 


2.297 


3.61 


.0003432 


0. 2662 


50 


7.85 


.000928 


.0000578 


1,596 


11 


3.203 


5.66 


.000555 


0.375 


31 


8.34 


.0009345 


.0000576 


12,730 


15 


4.105 


7. 78 


.001162 


0.466 


54 


15.55 


.00053 


.000415 


509 


12 


2.7065 


6.85 


.0002955 


0. 792 


58 


22.36 


.0006245 


.000485 


256 


♦10 


2.4085 


4.9 


.000495 


0.1137 


59 


22.0 


.0006246 


.000485 


2,045 


15 


3.311 


7.5 


.000848 


0. 1378 


61 


26.2 


.0006245 


.0003835 


324 


6 


2.511 


3.23 


.00114 


0.232 


63 


16.85 


.000524 


.0003315 


629 


8 


2.7985 


4.2 


.000566 


0.14 


64 


13.25 


.000533 


.0005355 


5,050 


13 


3.7035 


6.7 


.0009075 


0.178 



♦Obtained from test descriptions. 



TABLli III 

Calculations of the General Instability Stress 
by Ryder's Method 



Test No. 


A 5 


A f 


I, 


h 




K, 


K 










30 


.0490 


.0332 


.0006675 


.0000282 


.0205 


.330 


.1905 


.125 


.0144 


5,500 


.414 


31 


• 0474 


tt 


.000658 


n 


.0431 


.682 


.198 


.375 


.0433 


16,530 


.916 


32 


.0485 


n 


. 000664 


n 


.0826 


1.294 


.193 


.650 


.0750 


28,600 


1.177 


35 


.0646 


n 


.000815 


tt 


.0337 


.501 


.290 


.415 


.0477 


18,200 


1.595 


36 


.0640 


it 


. 000809 


n 


.0678 


1.012 


.292 


.615 


.0707 


27,000 


1.733 


37 


.0658 


ti 


.000823 


tt 


.1336 


1.965 


.284 


.934 


.1073 


41,000 


2.16 


38 


.0496 


tt 


.000675 


n 


.1625 


2.60 


.189 


.950 


.1092 


41,700 


1.55 


45 


.0604 


,0352 


.000795 


.0000412 


.0252 


.284 


.232 


.251 


.0296 


11,250 


.811 


46 


.0597 


n 


.00079 


n 


.0507 


.574 


.235 


.469 


.0551 


21,000 


1.325 


47 


.0598 


H 


. 00079 


tt 


.1014 


1.145 


.235 


.756 


.0888 


33,900 


1.433 


49 


.0856 


n 


.000938 


it 


.0427 


.400 


.328 


.478 


.0563 


21,500 


1.77 


50 


.0831 


n 


.000928 


M 


.0864 


.826 


.338 


.806 


.0948 


36,200 


2.54 


51 


.0822 


tt 


.000925 


H 


.1735 


1.670 


.342 


1.142 


.1346 


51,400 


2.16 


58 


.0456 


.0276 


.000625 


.000434 


.676 


.589 


.205 


1.538 


.180 


68,700 


1.665 


59 


.0456 


tt 


.000625 


tt 


1.352 


1.178 


.205 


1.934 


.226 


86,300 


lo475 


61 


.0456 


it 


.000625 


it 


.534 


1.937 


.205 


1.538 


.180 


68,700 


1.47 



TABLE IV 



Calculations of the General Instability Stress 
By Taylor 1 s liethod • 



Test 
86, 


*1 


h 


E 


l/tl 


l/t 2 




*• 






25 


. 0001419 


. 00000193 


.00383 


78.2 


274,0 


417.0 


.0228 


124,500 


7.75 


27 


.0001419 


. 00000771 


n 


78.2 


69.5 


it 


.0228 


147,000 


5.40 


30 


. 0000740 


. 00000193 


tf 


156,0 


274,0 


it 


,0164 


164,000 


12.80 


52 


• 0000740 


.00000771 


w 


156.0 


68.5 


tr 


,0164 


185,300 


7.3 


55 


.0000369 


« 00000193 


1! 


512.0 


274,0 


tr 


,0132 


186,000 


10.1 


37 


,0000369 


. 00000771 


tT 


312.0 


66.5 


it 


.0132 


209,000 


10*5 
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TABLE V 

Values of the Factor ^ and the blende rness Ratio \ 
for the Test Specimens 



Test No. 


i 


Pi 


via 


*5 




0 tl 
CO 


• liuo 


• UcOoO 


1 H77Q 

• iu/ / y 


a a K7 


Q*> KA 


CO 


TICK 

• lib o 


. ucb / ( 


• 1325 


0. 562 


56. 85 


c I 


• llbo 


• U2.oob 


1 COD 


0. b ( 0 


•7 C\ Af? 

39. 03 


CO 


« 110:7 


• Ulo / O 




a t£A 


JLO O* 4 


ou 


1 1 ft 
• 1 J.O o 






a ^ aa 


i£i. y 


u 1 


« llUJ 




• 1 1 c Crt 


A A77 


7fi QA 

/ o. yu 


6c 


• lib ( 


AOQ •? £J 


♦ 1345 


0. 570 


55. 28 


tA 


1 1 1 A 

• 111U 


• Ulo / o 


, Ubyo 


J. 2 54 


2 / y. 2 


3 5 


i t oa 
• 112U 


• 0233 5 


A«7 C>! 

07 54 


3.320 


175. 7 


36 


• 1122 


• Ocoii 


, 0926 


0.393 


116. 7 


37 


, 1117 


, 02836 


. 11192 


0,476 


79. 80 


38 


. 1164 


. 02801 


, 1593 


0.676 


39.39 


39 


. 1170 


• 02801 


1897 


0. 805 


27. 80 
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Fig. 1 



Note: All drawings are twice actual size 
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Fig. 10 
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Fig. 13 
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Fig. 14 




NACA Technical Note No. 909 



Fig. 15 
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Fig. 27 




NACA Technical Note No. 909 Fig. 29 




NACA Technical Note No ,909 




Technical note No. 909 




NACA Technical Note No* 909 Figs.. 43. 44 




Figure 44. Local instability failure of the channel section 
frames. 
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Figure 46. Typical panel instability type failure. 




Figure 48. Tension failure. 



